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Abstract 



Differential photoproduction cross sections are measured for events containing D*^ 
mesons. The data were taken with the HI detector at the ep collider HERA and correspond 
to an integrated luminosity of 51.1 pb _1 . The kinematic region covers small photon vir- 
tualities Q 2 < 0.01 GeV 2 and photon-proton centre-of-mass energies of 171 < W JP < 
256 GeV. The details of the heavy quark production process are further investigated in 
events with one or two jets in addition to the D*^ meson. Differential cross sections for 
Z)*+jet production are determined and the correlations between the Z)* ± meson and the 
jet(s) are studied. The results are compared with perturbative QCD predictions applying 
collinear- or /^-factorisation. 
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1 Introduction 



Charm photoproduction in ep collisions at HERA proceeds predominantly via photon-gluon 
fusion as shown in Fig. 1, where the quasi real photon (virtuality Q 2 ~ GeV 2 ) is emitted from 
the beam lepton. The charm quark mass provides a hard scale which allows perturbative QCD 
(pQCD) to be applied over the full phase space. Therefore, charm photoproduction is partic- 
ularly well suited to test perturbative calculations and the underlying theoretical approaches. 



a) Direct 7 b) Real Gluon Emission c) Excitation d) Hadron - like 




Figure 1: Diagrams for charm photoproduction: direct photon processes (a) and (b), and re- 
solved photon processes with (c) charm excitation and (d) the hadronic manifestation of the real 
photon. 



Previous measurements have focused on inclusive D* ± meson production [1], D*^ meson 
with associated dijet production [2-4] and heavy quark pair production using events with a 
D*^ meson and a muon [5]. Here D*^ photoproduction is considered, using a data sample five 
times larger than in the previous HI analysis [1]. Single and double differential cross sections 
for inclusive D* ± production are presented and compared with theoretical calculations using 
collinear-factorisation [6-9] or /c r factorisation [10-13]. Details of the heavy quark production 
process are investigated further by studying events either with an additional jet not containing 
the D* meson ("D* + other jet") or with two jets ("D* tagged dijet"). The jets are measured 
down to a transverse momentum of p t (jet) = 3 GeV, which extends the region explored in [2-4] 
to significantly smaller values. Whereas the D* always originates from a charm or anticharm 
quark produced in the hard subprocess, the non-D* -tagged jet can result from either the other 
heavy quark (Fig. la,b,d) or a light parton (e.g. a gluon, Fig. lb,c) coming from higher order 
processes. Measurements of correlations between the D* and the jet are performed which are 
sensitive to these higher order effects and to the longitudinal and transverse momenta of the 
partons entering the hard scattering process. 

The paper is organised as follows. In section 2 a brief description of the HI detector is 
given, followed by the details of the event selection and the reconstruction of the D* mesons. 
In section 4 the theoretical calculations are described. In section 5 the determination of the 
cross sections and the systematic uncertainties are presented and the measured cross sections 
are compared with theoretical calculations. 
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2 Detector Description 



The HI detector is described in detail in [14, 15] and only the components most relevant for this 
analysis are briefly mentioned here. The origin of the HI coordinate system is the nominal ep 
interaction point, the direction of the proton beam defining the positive ,2-axis (forward region). 
The transverse momenta are measured in the xy plane. 

Charged particles are measured in the Central Tracking Detector (CTD) which covers the 
range in pseudo-rapidity 1 1.74 > rj > —1.74. The CTD comprises two large cylindrical Central 
Jet Chambers (CJCs) arranged concentrically around the beam-line within a solenoidal mag- 
netic field of 1.15 T. Two additional drift chambers improve the ^-coordinate reconstruction. 
The track resolution is further improved using hit information from the central silicon track de- 
tector (CST) [16]. The CTD also provides trigger information which is based on measurements 
in the r-0 plane of the CJCs and the z-position of the interaction vertex obtained from a double 
layer of multiwire proportional chambers. 

The track detectors are surrounded in the forward and central directions by a fine grained 
Liquid Argon calorimeter (LAr, 3.4 > i] > —1.5) and in the backward region by a lead- 
scintillating fibre calorimeter (SpaCal, —1.4 > r] > —4.0) [17] both of which contain elec- 
tromagnetic and hadronic sections. For the present analysis the hadronic final state (HFS) is 
reconstructed from combined objects, built from calorimeter clusters and tracks, using an algo- 
rithm which ensures that no double counting of energy occurs. Compared to the case of clusters 
alone the use of combined objects improves the reconstruction of low momentum particles. 

A crystal Cerenkov calorimeter (electron tagger) located close to the beam pipe at z = 
—33.4 m detects positrons scattered through a very small angle (tt — 6 e i ) < 5 mrad and is used to 
trigger photoproduction events. The ep luminosity is determined from the QED bremsstrahlung 
(ep — > epq) event rate by detecting the radiated photon in another Cerenkov calorimeter located 
at z = — 103 m (photon detector). 

3 Event Selection and Reconstruction 

The data were recorded in e + p scattering at HERA in 1999 and 2000 and correspond to an 
integrated luminosity of C = 51.1 pb _1 . The energy of the positrons was E e = 27.6 GeV and 
that of the protons E p = 920 GeV. 

The events were triggered by requiring signals from the central drift chambers and the multi- 
wire proportional chambers in coincidence with a signal in the electron tagger. In addition, an 
on-line software filter selected events with candidates for charmed hadron decays by calculating 
invariant masses of track combinations. 

The analysis of photoproduction events is restricted to 0.29 < y < 0.65, where the inelas- 
ticity y = 1 — E e r/E e is calculated from the reconstructed positron energy E e > in the tagger. In 
this y-range the average tagger acceptance is almost 40%. This range corresponds to a photon- 
proton centre-of-mass energy 171 ^ W 7P ^256 GeV. For the small scattering angles considered 
the photon virtuality Q 2 is below 0.01 GeV 2 . 

'The pseudo-rapidity r\ corresponding to a polar angle 9 (measured with respect to the positive z-axis) is given 
by r? = - In tan(0/2). 
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Selection of 

D*+ -► D°TT+ -> K-7T+7T+ 


PtCK» > 0.3 GeV 
p t (7r s ) > 0.12 GeV 
20° < 6(K, 7T, tt s ) < 160° 
d£/dx(i i f, 7r, 7r s ) consistent with particle hypothesis 


m(Kn) — m D o < < 
Am = m( 


80MeV for Pt(D*) < 6.0 GeV 
100 MeV for 6.0 < Pt (D*) < 8.5 GeV 
140 MeV for pt{D ) > 8.5 GeV 

Ktttt,) - m(K7r) < 167.5 MeV 


Pt (K)+pt(n) > 2.2 GeV 
> 10 


Visible kinematic region 


inclusive D* 




Pt {D*) > 2.0 GeV 

\r](D*)\ < 1.5 
Q 2 < 0.01 GeV 2 
0.29 < y < 0.65 




D* + other jet 


D* tagged dijet 


|7?(jet)| < 1.5 
p t (jet) > 3 GeV 
-D* £ jet 


|77(jet)| < 1.5 
Pt(jeti) > 4 GeV 
Pt(jet 2 ) > 3 GeV 



Table 1: selection cuts and definition of the "visible" kinematic regions for which the inclu- 
sive D*, the D* + other jet and the D* tagged dijet cross sections are measured. 



3.1 D* Selection and Fit 

The -D*^ meson is detected via the decay channel 2 D* + — > D°nf — > K~n + irf, which has a 
branching ratio BTZ(D* — > A'tttts) = (2.57±0.06)% [18]. Here 7r s refers to the low momentum 
pion from the D* decay. In each event, tracks with opposite charges, fitted to the event vertex, 
are combined in pairs and both invariant masses m(K ± Tx T ) are calculated where one track is 
assigned the kaon mass and the other the pion mass. If the result is consistent with the nominal 
D° mass [18], any remaining track with a charge opposite to that of the kaon candidate is 
assigned the pion mass and combined with the D° candidate to form a D* candidate. The 
measured specific energy loss per path length dE/dx of each track has to be consistent with the 
respective particle hypothesis [19]. 

The D* candidate is accepted if it fulfils the selection cuts listed in Tab. 1. The restrictions 
in the transverse momentum pt(D*) and in the pseudorapidity rj(D*) ensure good detector ac- 
ceptance. To reduce the combinatorial background, cuts are applied to the scalar sum of Pt(K) 
and pt(n) and to the fraction of the transverse momentum carried by the D* with respect to the 
scalar sum of transverse energies of the full hadronic final state, excluding the forward region 
(9 < 10°). 

In Fig. 2a the distribution of the mass difference Am = m(Kirn s ) — m(Kir) of the fi- 
nal D* candidates is shown. A clear peak is observed around the nominal value of Am = 

2 Charge conjugate states are implicitly implied. 
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(145.421 ± 0.010) MeV [18]. The number of reconstructed D* mesons N(D*) is extracted 
from a likelihood fit to the Am-distribution with a function which is a superposition of a 
Gaussian for the signal and a phase space threshold function with a quadratic correction term, 
F b (Am) = u n (Am — m^)"" 2 • (1 — w s (Am) 2 ), for the background. The fit yields a signal of 
1 166 ± 82 D* mesons. Separate fits are performed in each bin of the cross section measurement, 
in which the mean and the width of the Gaussian, the background parameters u s and, in cases 
where statistics are very low, u e are fixed to the values obtained from the fit to the inclusive D* 
sample. The background normalisation u n is a free parameter. 

3.2 Jet Selection 

In the inclusive D* sample jets are defined by the inclusive /c t -algorithm [20] in the ^-recom- 
bination and Ai?-distance scheme (with R = 1, [21]). The jet algorithm is applied in the 
laboratory frame to all reconstructed particles of the HFS. To prevent the decay particles of the 
D* candidate from being ascribed to different jets, the D* candidate is used as a single particle 
in the jet algorithm, replacing its decay products. In events which contain more than one D* 
candidate, the jet algorithm is run separately for each candidate. The transverse momentum 
is required to satisfy p*(jet) > 3 GeV. For jets originating from charm the requirement of 
Pt(jet) > 3 GeV matches approximately the requirement of Pt(D*) > 2 GeV, since the charm 
fragmentation function to a D* peaks where the D* meson takes ~ 70% of the charm quark 
momentum [22]. 

To ensure a good jet reconstruction at these low values of p t , jets are restricted to the central 
detector region |?y(jet)| < 1.5 where precise track information is available. The jet transverse 
momentum resolution is around 30% over the whole momentum range considered here. 

The jet with the highest^ (jet) not containing the D* meson is considered together with the 
D* meson in the ll D* + other jet" analysis. The Am-distribution of the D* candidates of this 
sample is shown in Fig. 2b. A signal of 592 ± 57 D* mesons is found. In about 10% of the D* 
events a second (non-/J*-)jet is observed. 

In addition to the D* + other jet sample a ll D* tagged dijet" sample is selected. At least 
two jets are required with £>t(jet) > 4 GeV for the highest p t jet and p*(jet) > 3 GeV for the 
second jet (see Tab. 1), irrespective of whether the D* meson is attributed to either of the jets. 
A signal of 496 ± 53 D* tagged dijet events is observed. The fraction of events where the D* 
meson is contained in neither of the two jets is negligible. 

4 QCD Calculations 

Figure 1 shows examples of diagrams for charm photoproduction in leading order 0(a s ) and 
next-to-leading order 0(a 2 s ). In direct-photon processes (Figs. la,b), the photon emitted from 
the beam lepton enters the hard process directly. In resolved-photon processes (Figs. lc,d), the 
photon acts as a source of incoming partons, one of which takes part in the hard interaction. 
The distinction between these two classes of processes depends on the factorisation scheme and 
the order in which the calculation is performed. 
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The data presented in this analysis are compared with leading order (LO) calculations sup- 
plemented by parton showers as well as with next- to-leading order (NLO) calculations. The 
calculations are performed using either the collinear factorisation or the /^-factorisation ap- 
proach. The collinear factorisation makes use of the DGLAP [6-9] evolution equations. In this 
approach transverse momenta obtained through the initial state QCD evolution are neglected 
and all the transverse momenta are generated in the hard scattering process, i. e. the incoming 
partons are collinear with the proton, resulting at lowest order in a back-to-back configuration 
of the heavy quark pair in the transverse plane. Effects from the finite transverse momentum 
of the gluons enter only at the NLO level (for example in a process like — > ccg as shown 
in Fig. lb). In the /^-factorisation ansatz the transverse momentum (k t ) of incoming gluons are 
already included at leading order both in the /^-dependent off-shell matrix element and the /in- 
dependent unintegrated gluon density [23]. Therefore, higher order corrections, i. e. hard parton 
emissions, are partially considered. 

Heavy quark production is calculated either in the massive scheme [24], where heavy quarks 
are produced only perturbatively via boson gluon fusion, or in the massless scheme [25], where 
heavy quarks are treated as massless partons. These two schemes should be appropriate in dif- 
ferent regions of phase space [26]: the massive scheme should be reliable when the transverse 
momentum p t of the heavy quarks is of similar size compared to m c , whereas the massless 
scheme is expected to be valid for p t ^> m c . Recently, new calculations combining and match- 
ing the two approaches in the photoproduction regime have become available [27, 28]. 

The uncertainties on the calculations are estimated by varying the charm quark mass, the 
renormalisation scale and where possible the factorisation scales. The uncertainty in each bin 
is obtained by taking the maximal deviations from the central value resulting from the separate 
variations. The detailed settings are summarised in Tab. 2. All predictions (except where scale 
dependent fragmentation functions are used) are based on a D* fragmentation ratio of f(c — > 
D*) = 0.235 [35]. The beauty contribution is included in all predictions and amounts to a few 
percent. 

PYTHIA In Pythia [36] three different processes are generated separately using leading or- 
der matrix-elements: direct photon-gluon fusion (Fig. la), resolved photon processes in which 
a charm quark (charm excitation, Fig. lc) or a gluon (Fig. Id) from the photon enters the hard 
scattering. In the excitation processes the charm quark is treated as a massless parton, whereas 
in the other processes the charm mass is accounted for in all steps of the calculation. Higher 
order contributions are simulated with leading log parton showers in the collinear approach. 
The Lund string fragmentation model [37] is used for the simulation of the hadronisation pro- 
cess. For the longitudinal fragmentation of the charm quark into the D* meson the Peterson 
parametrisation [38] is used. No uncertainties are calculated. 

CASCADE The CASCADE [39-41] program is used for leading order calculations in the 
/^-factorisation approach. In the ^g* — > cc matrix element, which takes the charm mass into 
account, the incoming gluon is treated off mass-shell and can have a finite transverse momen- 
tum. Higher order QCD corrections are simulated with initial state parton showers applying 
the CCFM evolution [42-45] equations with an unintegrated parton density function (uPDF) 
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PYTHIA CASCADE FMNR ZMVFNS GMVFNS 



Version 6.224 1.2010 



Proton PDF CTEQ6L [29] A0(±) [30] CTEQ6M [29] CTEQ6M CTEQ6M 

Photon PDF GRV-GLO[31] - GRV-GHO[31] GRV-G HO GRV-G HO 



2 1 

Renorm. scale m t 1 

0.5 J 


\-m' t 1 
1 0.5 




i 'I 

> ■ m t ,D* 1 , 

0.6 J 


> ■ m t ,D* 


1 

Factor, scale m t y/s + Q\ 2 > 

4j 


1/1.5 ] 


0.6 | 

> ■ m t ,D* 1 
2 J 


> ■ m t ,D* 



m c [GeV] 1.5 1.5tg-j 1.5±^;| ~ 1.5 1.5 



Fragmentation e pet = 0.04 e pet = 0.04 e pet = 0.035 [32] BKK O [33] KKSS [34] 

Table 2: Parameters used in the pQCD calculations where m c denotes the charm quark mass. 
The "transverse mass" variables are defined as m 2 = ml + {p 2 c + pf 5 )/2, m 2 = 4m 2 c + p 2 c 
and m 2 Df . = m 2 c + p 2 D * where p tjC (p tjS ) is the transverse momentum of the charm quark 
(antiquark) andp tt o* is the transverse momentum of the D* meson. The squared invariant mass 
and the transverse momentum squared of the cc pair are denoted by s and Q\, respectively. e pet 
is the Peterson fragmentation parameter. If a parameter is varied to determine the uncertainty 
of the prediction, the central line gives the value used for the main prediction. 



including angular ordering constraints for the emitted partons. The uPDF has been obtained 
from an analysis of the inclusive structure function F 2 in the CCFM approach [30]. For the 
variation of the renormalisation scale, dedicated unintegrated gluon density parameterisations 
have been used [30]. The factorisation scale cannot be varied. The final state radiation off the 
heavy quarks and the fragmentation is performed with Pythia. 

FMNR The FMNR program [24, 46] implements a massive scheme next-to-leading order 
calculation (0(a 2 )) in the collinear factorisation approach. It provides weighted parton level 
events with two or three outgoing partons, i. e. a cc quark pair and possibly one additional light 
parton. The fragmentation of the charm quarks into D* mesons is treated by a downscaling 
of their three-momenta according to the Peterson fragmentation function in a frame where the 
quark-antiquark pair is back-to-back. 

ZMVFNS The zero-mass variable-flavour-number scheme (ZMVFNS) [25, 47] provides a 
next-to-leading order calculation (0(a 2 )) for D* + other jet cross sections in the collinear ap- 
proach, neglecting the charm mass. The transition from the charm quark to the D* meson is 
treated using the scale dependent fragmentation functions determined in [33]. In the determi- 
nation of the uncertainty the initial and final state factorisation scales are varied simultaneously 
(with a lowest value of m t: D* and 1.5 m t: D*, respectively, see Tab. 2). 

GMVFNS The general-mass variable-flavour-number scheme (GMVFNS) combines the mass- 
less with the massive scheme [27,28]. Scale dependent fragmentation functions as determined 
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in [34] are used. The calculation is only available for inclusive D* production. The initial and 
final state factorisation scales are varied separately in the uncertainty determination. 

For the NLO calculations which are compared to the D* + other jet and D* tagged dijet 
measurements (FMNR and ZMVFNS), additional corrections for hadronisation effects (tran- 
sition form partons to jets) are applied. These hadronisation corrections are calculated using 
PYTHIA. In Pythia parton level jets are constructed from the generated quarks and gluons af- 
ter the parton showers have been simulated. The ratio of hadron and parton level cross sections 
in each bin is applied to the NLO calculations as a hadronisation correction factor. 



5 Results 

5.1 Cross Section Measurement 

The bin averaged visible differential cross section with respect to a variable Y (with bin width 
AY) is calculated according to 



where a correction r = 0.035 [1] is applied to account for reflections from other D° decays 
within the D° mass window. B1Z(D* — > Kirir s ) is the branching ratio of the analysed D* decay 
chain and C is the integrated luminosity. The correction factor e takes into account detector 
acceptances, trigger and reconstruction efficiencies and migrations between bins. The cross 
section is defined as the sum of the D* + and D*~ cross sections and includes D* mesons from 
6-quark decays. 

The average acceptance of the electron tagger is calculated as a convolution of the pre- 
dicted y - distribution with the y - dependent tagger acceptance which is determined as in [48]. 
The reconstruction and trigger efficiencies as well as the acceptance of the selection criteria 
are determined using a GEANT 3.15 [49] based simulation of the detector response to events 
generated with Pythia. The efficiency of the particle identification using dE/dx and of the 
software filter are deduced from the data. 

The following sources of systematic uncertainties (summarised and quantified in Tab. 3 for 
the integrated cross section) have been studied: 

• The simulation of the trigger signals from the CJCs and the multi-wire proportional cham- 
bers have been verified (using data) to within an uncertainty of 4%. Including the ef- 
ficiency of the software filter this amounts to a 4.5% uncertainty on the cross section 
measurement. 

• The positron beam parameters and the absolute energy scale of the electron tagger have 
been varied within their uncertainties. The resulting average uncertainty on the tagger 
acceptance is 5.8%, depending on W 1P . 

• The uncertainty on the track reconstruction efficiency amounts to 6% per D* meson. 



da vis (ep -> e> 'D^j+jeQX) 
dY 



AY ■ BTZ(D* -> Kmi s ) -C-e 



iV(D* ± (+jet)) • (1 - r) 



(1) 
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sources 


inclusive D* 


D* with jets 


[%] 


[%] 


trigger efficiency 


4.5 


electron tagger acceptance 


5.8 


track reconstruction 


6 


particle identification 


2 


&Jl(D* -> Kitit s ) 


2.5 


reflections 


1.5 


signal extraction 


3 


model dependence 


1.3 


luminosity C 


1.5 


HFS objects energy scale 




2.8 


in total 


10.7 


11.1 



Table 3: Systematic uncertainties of the integrated cross section measurements. 



• The uncertainty on the efficiency of the particle identification using dE/dx is estimated to 
be 2% per D* meson [19]. 

• The uncertainty on the D* branching ratio is 2.5% [18]. 

• The uncertainty on the correction for reflections from other D° decays is 1.5% [1]. 

• The uncertainty on the extraction of the D* signal from the Am distributions has been 
determined in [19] and amounts to 3%. 

• The model dependence of the correction factors is estimated by using CASCADE instead 
of PYTHIA. Half of the resulting deviation is taken as the uncertainty on the signal ex- 
traction procedure, i.e. 1.3% for the total inclusive D* sample and up to 17% for the 
differential distributions. 

• The luminosity is determined with a precision of 1.5%. 

• The uncertainty on the energy of the HFS objects leads to a systematic uncertainty on the 
D*-jet cross sections of 2.8%. 

The total systematic uncertainty is obtained by adding each uncertainty in quadrature and 
amounts to 1 1% for the integrated cross sections and up to 20% for differential distributions. 



5.2 Inclusive D* Cross Sections 



The integrated D* photoproduction cross section is measured to be: 

a vis (ep -> e'D** X) = 6.45 ± 0.46 (stat.) ± 0.69 (sys.) nb, 

in the visible range given by 0.29 < y < 0.65, Q 2 < 0.01 GeV 2 , p t (D*) > 2 GeV and 
\r)(D*)\ < 1.5. In Tab. 4 the predictions from QCD calculations are listed and compared with 
this measurement. The central values from FMNR and CASCADE are slightly lower than the 
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o 'vis [nb] 


inclusive D* 


D* + other jet 


D* tagged dijet 


Data 


6.45 ± 0.46 ± 0.69 


3.01 ± 0.29 ± 0.33 


2.32 ± 0.25 ± 0.26 


FMNR 




(2.65l°;i) 


(2.44l°f) 


® had. corr. 




o qpr+0.69 
z -°°-0.37 


2.09±g;S 


ZMVFNS 


- 


(3.051H?) 


- 


® had. corr. 




9 71 +0.55 
z - ' 1 -0.42 




GMVFNS 


8-2^ 






Pythia 


8.9 


3.8 


2.8 


Cascade 


r OO+0.54 
°-°°-0.62 


3.08lg;i| 


2.481™ 



Table 4: Integrated cross section in the visible range for inclusive D*, D* + other jet and D* 
tagged dijet photoproduction. For the cross sections for processes involving jets the result of 
the NLO calculations before the correction for hadronisation effects is given in brackets. 



measured result, whereas those of Pythia and GMVFNS are higher. 

The measured bin averaged differential cross sections are shown in Figs. 3-5 and given in 
Tabs. 5-6. The predictions show large uncertainties compared with those on the data. These 
uncertainties partially cancel in the comparison of the shapes of the distributions. The ratio 

1 da calc 

calc Ay 

T3 _ a vis 111 rrj\ 
1 rfqdata W 
data d y 

is therefore also presented, where Y denotes any measured variable. Note that here each differ- 
ential cross section is normalised to its own visible cross section (inclusive D*, D* + other jet 
or D* tagged dijet). 

In Fig. 3a and b the cross section is shown differentially in the transverse momentum of 
the D*. The cross section falls steeply with increasing p t (D*) as predicted by all calculations. 
CASCADE and FMNR predict a distribution which falls less steeply at large p t (D*) than is vis- 
ible in the data. The PYTHIA and GMVFNS calculations describe the slope well. However, the 
theoretical uncertainty of the GMVFNS calculations is large. In Fig. 3b the small contribution 
coming from 6-decays as determined from FMNR is shown separately. 

In Fig. 3c and d the differential cross section is shown as a function of the pseudorapidity 
r)(D*). The cross section decreases with increasing 77 (the forward direction 77 > 0). All 
calculations predict a similar shape, differing from that in the data, which shows a larger relative 
contribution in the forward direction. 

In Fig. 4 the differential cross section is shown as a function of r)(D*) for three bins in 
p t (D*). As for the inclusive rj(D*) distribution none of the calculations is able to describe the 
data completely. 

The inelasticity z(D*), defined by z(D*) = P ■ p(D*) / (P ■ q) withP,p(D*) andgbeingthe 
four-momenta of the incoming proton, the D* meson and the exchanged photon, is a measure of 
the fraction of photon energy transferred to the D* meson in the proton rest frame. This quantity 
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is sensitive to the production mechanism and to the c — > D* fragmentation function. It is recon- 
structed as z(D*) — (E — p z ) D * I (2yE e ). The differential cross section as a function of z(D*) 
is compared to the predictions in Fig. 5a and b. For z(D*) > 0.2 the CASCADE and FMNR 
predictions agree fairly well with the data. However, in both cases the predictions (central val- 
ues) underestimate the measured cross section in the region z(D*) < 0.2. A similar finding is 
reported for a measurement of D* production in deep inelastic scattering (Q 2 > 1 GeV 2 ) [50]. 
PYTHIA can account for the large cross section at small z(D*), but overestimates the cross sec- 
tion at medium z(D*). The central value of the GMVFNS calculation is slightly higher than the 
measurement, but the uncertainties cover the whole range. None of the calculations describes 
the overall shape of the data. 

The cross section as a function of W rp is shown in Fig. 5c and d. It falls weakly with 
increasing W 1P . The behaviour observed in the data is reproduced by all calculations within the 
uncertainties. 

In summary, significant differences are observed in the shape of most of the distributions 
between data and all central theoretical predictions. A good agreement is only found for the 
W 7P distribution for which each bin receives contributions throughout the available phase space 
of the produced D* ± meson. Here all calculations predict very similar shapes. For those ob- 
servables which are especially sensitive to the phase space distribution of the outgoing charm 
quark, i.e. p t (D*), r/(D*) and z(D*), large deviations in shape between data and theory are 
observed. It is interesting to observe that the cross sections predicted by the calculations which 
treat the heavy quark mass explicitly (CASCADE and FMNR) are similar. 

5.3 D*+ other Jet and D* tagged Dijet Cross Sections 

A more detailed investigation of the heavy quark production process is performed by analysing 
events with a D*^ meson with either a jet not containing the D* (D* + other jet ) or with 
two jets (D* tagged dijet ). In this way one can tag a second outgoing parton from the hard 
interaction in addition to the (anti-)charm quark. The requirements for the D* meson and the 
photoproduction selection are the same as in the inclusive case. In the D*+ other jet analysis 
the jet is required to have a transverse momentum p t (jet) > 3 GeV in the range of \r) (jet) | < 1.5 
in the laboratory frame. The integrated (D* + other jet) cross section in the visible range given 
in Tab. 1 is measured to be 

Vvis{ep -»• e' J D* ± + other jet X) = 3.01 ± 0.29 (stat.) ± 0.33 (sys.) nb. 

The integrated D* tagged dijet cross section (pt(jet) > 4(3) GeV, |?y(jet)| < 1.5) is measured 
to be (in the visible range given in Tab. 1) 

Vvis(ep -> e'D^+ljet X) = 2.32 ± 0.25 (stat.) ± 0.26 (sys.) nb. 

In Tab. 4 both values are compared with the predictions from the QCD calculations. All pre- 
dictions agree with the measurements within the quoted calculated uncertainties, including 
ZMVFNS for the D* + other jet sample. 

The bin averaged differential cross section for the D* + other jet and D* tagged dijet are 
listed in Tabs. 7-8 and shown in Figs. 6-9. The cross sections for the D* + other jet selection 
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are shown in Fig. 6 as functions of Pt(D*) and p t (jet) together with theoretical predictions. The 
calculations overestimate the cross section at large p t (D*), while the cross section as a function 
of p t (jet) is well described by the NLO calculations. CASCADE overestimates the cross sections 
at large p t (jet). 

Cross sections as a function of r)(D*) and r?(jet) are shown in Fig. 7a-d. They differ notice- 
ably: The r](D*) distribution falls steeply with increasing values of r\ (similar to the inclusive 
case), whereas ^(jet) is almost flat. For direct photon-gluon fusion processes (^g — > cc), as 
shown separately for Pythia, the shapes of the cross sections as a function of r](D*) and ^(jet) 
are found to be similar (dotted lines in Fig. 7a and c), indicating that the difference is not caused 
by the slightly different kinematic cuts for the D* and the jet. The observed difference in shape 
between the r/(D*) and r?(jet) distributions indicates the presence of hard non-charm partons in 
the forward region. In fact, the dominant mechanism for non-charm jets, as predicted by the 
calculations, is hard gluon radiation off the gluon from the proton, either calculated as an NLO 
correction (Fig. lb), as a (LO) charm excitation (Fig. lc) or by using an unintegrated gluon 
density. All calculations include this diagram and can describe the observed shapes of 7](D*) 
and rj (jet) in general. 

The D* + other jet measurement allows the investigation of correlations between the D* 
and the jet. The cross section as a function of i](D*) -r?(jet) is shown in Fig. 7e and f. The 
distribution is not symmetric but the jet is on average more forward than the D*. The cross 
section is reasonably well described in shape and magnitude by all QCD calculations. 

In collinear factorisation at LO a resolved photon process is characterised by x 7 < 1, where 
Xj is the fraction of the longitudinal photon momentum entering the hard scattering process. In 
the D* tagged dijet sample, x 7 is experimentally approximated by 

obs _ Z^i=l Kf^ Pzji 

The sum in the numerator includes the particles in the two selected jets, whereas the sum in 
the denominator contains all reconstructed particles of the hadronic final state. At NLO x° hs 
is sensitive to 0(a 2 s ) contributions. In the /^-factorisation approach the rr° bs observable is 
sensitive to the contribution from gluon emission in the initial state. 

The D* tagged dijet cross section as a function of x° bs is shown in Fig. 8. The large cross 
section at small x° shows that processes beyond direct photon-gluon fusion are needed to 
describe the data in the collinear approach (as can be seen by comparison with the prediction 
from Pythia (dir.) in Fig. 8a.). Both Pythia and Cascade give a poor description of the 
x° hs distribution. All predictions underestimate the region of low x° hs < 0.6 which can be 
clearly seen in the normalised shape R. The hadronisation corrections applied to the FMNR 
calculation are large. 

The correlation in the transverse plane is experimentally accessed by the difference in the 
azimuthal angle A0(D* jet) between the D* and the jet in the D* + other jet sample. The cross 
section as a function of A(f)(D*, jet) is shown in Fig. 9. Only ~ 25% of the measured cross 
section originates from a back-to-back configuration in the transverse plane with A(f)(D*, jet) > 
170°. Such a configuration is expected for the process — > cc in the collinear approximation. 
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The large fraction of events where the D* and the jet are not back-to-back can only be described 
by models which include significant contributions from higher order QCD radiation. 

In Pythia higher order contributions are simulated by leading log parton showers and the 
charm excitation process. In CASCADE higher order contributions lead to a significant k t of the 
gluon entering the hard subprocess as simulated by the unintegrated gluon density. The large 
deviations from the back-to-back topology are well described by both PYTHIA and CASCADE. 
The back-to-back configuration is overestimated by Pythia, while with the given parameter- 
isation of the unintegrated gluon density CASCADE tends to underestimate the back-to-back 
region and overestimates the small A(f)(D*, jet) configuration which corresponds to large k t . It 
is interesting to note that CASCADE also predicts significantly harder p t spectra than that ob- 
served in the data for the inclusive D* cross section in Fig. 3a as well as for the D* and jet cross 
sections in Fig. 6a and c for the D* + other jet sample. These quantities are also sensitive to the 
k t of the gluon which suggests that the current unintegrated gluon density is too hard in k t . In 
the NLO calculations A0(D* jet) 7^ 180° comes entirely from real gluon emission (7$ — > ccg) 
and from processes initiated by light quarks from the proton (75 — > ccq). The NLO calculations, 
although 0(a 2 s ), include only the lowest order contribution to this region. The NLO calcula- 
tions are in reasonable agreement with the measurement in the region of A<fi(D*, jet) > 120°, 
but show large discrepancies for A(j)(D*, jet) < 120°, suggesting the presence of higher order 
contributions. A similar behaviour was found in the measurement of D* dijet data by the ZEUS 
collaboration [4] where significantly higher transverse momenta of the jets were required. 



6 Conclusions 

The photoproduction of D* mesons is investigated with the HI detector at HERA using a data 
sample five times larger than in a previous publication. Differential cross sections are deter- 
mined for events with a D* meson (inclusive D*) and for events with a D* meson and jets 
(D* + other jet and D* tagged dijet). Jets are selected with transverse momenta down to 3 GeV. 
The results are compared with QCD calculations based on different approaches in leading and 
next-to-leading order pQCD employing either collinear or /c r factorisation. 

The precision of the inclusive D* cross section measurements presented here is much higher 
than the accuracy of the NLO calculations. The comparison of normalised cross sections, for 
which these theoretical uncertainties are significantly reduced, reveals sizable differences be- 
tween data and theoretical predictions for variables sensitive to the phase space distribution of 
the outgoing charm quark. However, a good agreement is found for the W 1P distribution for 
which each bin receives contributions throughout the available phase space of the produced 
meson. 

In the D* + other jet analysis a jet is required in addition which does not contain the D* 
meson. In general the predictions agree better with the measurements than is the case for the in- 
clusive D* analysis. The cross section as a function of r?(jet) is almost flat in contrast to the cross 
section as a function of r](D*) which falls towards increasing r)(D*). This indicates the presence 
of jets originating from non-charmed partons, most likely from gluons. All calculations are able 
to reproduce this feature of the data. Correlations between the D* and the jet in the transverse 
plane are investigated by measurements of the difference in azimuthal angle A(f)(D*, jet). A 
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large fraction of the produced _D*+jet combinations deviates from a back-to-back configuration, 
indicating the importance of higher order contributions. The large tails of the A(p(D*, jet) dis- 
tribution can be reproduced reasonably well by LO calculations in the collinear factorisation 
ansatz which include parton showers and charm excitation processes and by using unintegrated 
gluon densities in the /c r factorisation ansatz. However, there are differences in shape observed 
between the data and these predictions. Especially for CASCADE the A<j)(D*, jet) distribution, 
together with the distributions in Pt(D*) and p t (jet), suggests that the unintegrated gluon den- 
sity used for this analysis overestimates the high k t region. The available NLO calculations 
underestimate significantly the observed cross sections in the region A(f)(D*, jet) < 120°. 

In the D* tagged dijet analysis two jets are required in addition to the D* meson and the 
observable a;° bs is studied. All calculations underestimate the cross section in the region of low 
£° bs < 0.6 where resolved photon processes or other higher order contributions are expected to 
be enhanced. 
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Inclusive D* Cross Sections 


p t (D*) range 


da/dp t (D*) 


stat. 


sys. 


[GeV] 


[nb/GeV] 




2 


2 5 


3 37 


63 


0.36 


2 5 


3 


3 01 


37 


0.32 


3 


3 5 


2 22 


23 


0.24 


3.5 


4.25 


0.97 


0.11 


0.10 


4.25 


5.0 


0.577 


0.071 


0.063 


5.0 


6.0 


0.250 


0.038 


0.027 


6.0 


8.5 


0.063 


0.013 


0.0068 


8.5 


12.0 


0.0133 


0.0050 


0.0016 


n( D*) 


range 


da/dn(D*) 


stat. 


sys. 












-1 5 


-1 


3 54 


30 


0.38 


-1.0 


-0.5 


2.93 


0.25 


0.31 


-0.5 


0.0 


1.71 


0.26 


0.19 


0.0 


0.5 


2.16 


0.29 


0.23 


0.5 


1.0 


1.58 


0.31 


0.18 


1.0 


1.5 


0.74 


0.41 


0.08 


z(D*) 


range 


da/dz(D*) 


stat. 


sys. 








[nb] 







0.1 


10.8 


2.5 


1.2 


0.1 


0.2 


17.8 


1.9 


2.0 


0.2 


0.3 


9.9 


1.2 


1.1 


0.3 


0.45 


8.66 


0.76 


0.93 


0.45 


0.6 


6.00 


0.58 


0.64 


0.6 


0.8 


1.48 


0.24 


0.19 


W 1P range 


da/dW lp 


stat. 


sys. 


[GeV] 


[nb/GeV] 




172 


192 


0.112 


0.012 


0.017 


192 


212 


0.0784 


0.0077 


0.0077 


212 


232 


0.0601 


0.0071 


0.0058 


232 


256 


0.0548 


0.0078 


0.0057 



Table 5: Bin averaged differential cross sections for inclusive D* production in bins ofp t (D*), 
rj{D*), z(D*) and W 1P with their statistical and systematical uncertainties. 



19 



Inclusive D* Cross Sections 


2.0 < p t (D*) < 3.0 GeV 


i](D*) range 


da/di](D*) [nb] stat. sys. 


-1.5 -0.9 
-0.9 0.2 
0.2 1.5 


2.14 0.25 0.23 
1.14 0.17 0.12 
0.63 0.20 0.07 


3.0 < pt(D*) < 4.5 GeV 


r)(D*) range 


da/dr](D*) [nb] stat. sys. 


-1.5 -0.9 
-0.9 0.2 
0.2 1.5 


1.02 0.11 0.11 
0.714 0.073 0.077 
0.495 0.082 0.057 


4.5 < p t {D*) < 8.5 GeV 


r](D*) range 


da/dr](D*) [nb] stat. sys. 


-1.5 -0.9 
-0.9 0.2 
0.2 1.5 


0.263 0.043 0.028 
0.240 0.037 0.027 
0.190 0.034 0.020 



Table 6: Bin averaged differential cross sections for inclusive D* production in bins ofr)(D*) 
for three ranges in Pt(D*) with their statistical and systematical uncertainties. 
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D*+ other Jet Cross Sections 




p t (D*) range 


da/d Pt (D*) 


stat. 


sys. 


J had 


[GeV] 


[nb/GeV] 






2.0 3.5 


1.35 


0.16 


0.15 


0.87 


1 C C A 
J.J J.U 


0.407 


0.051 


A A/1 A 

U.U44 


0.92 


SO 7 S 

J.U I.J 


0.124 


0.017 


014 


O Qf> 


7 S 110 
/ . J ll.u 


0.0116 


0.0048 


O 001 'X 


0.97 


Ptijet.) range 


da/dp t (}et) 


stat. 


sys. 


J had 


[GeV] 


[nb/GeV] 






3.0 5.0 


0.870 


0.094 


0.095 


A A 1 

0.91 


5.0 7.5 


0.334 


0.049 


0.037 


0.85 


/.j 11. U 


A AO A 


A AO 1 


A A1 A 


0.85 


110 16 
1 l.U 1 u.u 


A AOA1 

U.UzUl 


A AAT 1 
U.UU/ 1 


O 0099 


Q1 


r)(D *) range 


da/dr](D*) 


stat. 


sys. 


f 

J had 




[nb] 








-1.5 -1.0 


1.46 


0.18 


0.16 


0.86 


-1.0 -0.5 


1.46 


0.17 


0.16 


A OO 

0.89 


-0.5 0.1 


0.99 


0.16 


0.11 


0.90 


A 1 AO 
U. 1 U.O 


1.05 


0.18 


A 1 O 


0.91 


8 1 S 
u.o i.J 


0.33 


0.22 


O 04 


vj.yvj 


r? (jet) range 


da/drj (jet) 


stat. 


sys. 


J had 




[nb] 








-1.5 -1.0 


0.84 


0.16 


0.092 


0.66 


-1.0 -0.5 


1.15 


0.19 


0.13 


0.90 


-0.5 0.1 


1.15 


0.18 


0.13 


0.94 


0.1 0.8 


0.84 


0.17 


0.093 


0.99 


AO 1 C 

0.8 1.5 


0.99 


0.17 


A 1 1 

0.1 1 


0.94 


r](D*) -r?(jet) range 




stat. 


sys. 


J had 




[nb] 








-2.8 -1.9 


0.318 


0.065 


0.037 


1.01 


-1.9 -1.0 


0.66 


0.10 


0.07 


0.98 


-1.0 0.0 


1.01 


0.13 


0.11 


0.86 


0.0 1.0 


0.63 


0.13 


0.07 


0.81 


1.0 1.9 


0.44 


0.11 


0.05 


0.88 


1.9 2.8 


0.084 


0.071 


0.012 


0.89 



Table 7: Bin averaged differential cross sections for the D* + other jet sample in bins ofpt(D*), 
Ptijet), r](D*), rjijet) and f](D*) - r\{jet) with their statistical and systematical uncertainties. 
The last column shows the hadronisation correction factors applied to the NLO calculations. 
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D* tagged Dijet Cross Sections 




x° bs range 


da/dx° bs 


stat. 

[nb] 


sys. 


fhad 


0.0 0.3 


0.60 


0.39 


0.12 


0.69 


0.3 0.6 


2.23 


0.45 


0.36 


0.88 


0.6 0.8 


2.12 


0.51 


0.27 


1.40 


0.8 1.0 


4.81 


0.44 


0.53 


0.68 


D*+ other Jet Cross Sections 




A(f>(D*, jet) range 


da/dA(f)(D*,jet) 


stat. 


sys. 


fhad 


n 


[nb/°] 






86 


0.00167 


0.00086 


0.00022 


0.62 


86 114 


0.0108 


0.0029 


0.0012 


0.91 


114 138 


0.0156 


0.0040 


0.0018 


0.91 


138 154 


0.034 


0.0062 


0.0039 


0.94 


154 170 


0.0561 


0.0077 


0.0061 


0.90 


170 180 


0.073 


0.010 


0.008 


0.91 


154 180 


0.0625 


0.0061 


0.0068 


0.91 



Table 8: Bin averaged differential cross sections for the D* tagged dijet sample in bins ofx° 
and the D* + other jet cross section as a function of A(f)(D*,jet) with their statistical and 
systematical uncertainties. The last column shows the hadronisation correction factors applied 
to the NLO calculations. 
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Figure 2: D* signal in the distribution of the mass difference Am = m(KiriT s ) — m(K7i) of 
the inclusive D* (a) and the D* + other jet (b) selection. The solid lines represent the fits to 
determine the number of D* mesons in the signal and the dashed lines indicate the resulting 
background parametrisations. 
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Figure 3: Inclusive D* cross sections as a function ofp t (D*) (a,b) and r)(D*) (c,d) compared 
with the predictions of PYTHIA and CASCADE on the left and of the next-to-leading order 
calculations FMNR and GMVFNS on the right. For FMNR the beauty contribution is shown 
separately for p t (D*). Here and in the following figures the inner error bars indicate the statis- 
tical errors and the outer error bars show the statistical and systematic uncertainties added in 
quadrature. The normalised ratio R ( Eq. 2 ) is also shown. 
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Figure 4: Inclusive D* cross sections as a function ofr](D*)for three bins ofp t (D*) compared 
with the predictions of PYTHIA and CASCADE on the left and of the next-to-leading order 
calculations FMNR and GMVFNS on the right. 
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Figure 5: Inclusive D* cross sections as a function of the inelasticity z(D*) (a,b) and the 
photon-proton centre-of-mass energy W 1P (c,d) compared with the predictions of PYTHIA and 
CASCADE on the left and of the next-to-leading order calculations FMNR and GMVFNS on 
the right. 
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Figure 6: £>* + other jet cross sections as a function of the transverse momenta of the D* and 
the jet compared with the predictions of PYTHIA and CASCADE on the left and of the next-to- 
leading order calculations FMNR and ZMVFNS on the right. Here and in the following figures 
the central FMNR prediction is shown before and after applying the hadronisation corrections 
for the jet. 
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Figure 7: D* + other jet cross sections as a function of the i](D*) and rj(jet) and of their 
difference, compared with the predictions of PYTHIA and CASCADE (left) and of the NLO 
calculations FMNR and ZMVFNS (right). Here and in the following figures the direct photon 
contribution of PYTHIA is shown separately and labelled "Pythia 6.2 (dir.)". 
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Figure 8: D* tagged dijet cross sections as a function of x°^ s compared with the predictions of 
PYTHIA and CASCADE on the left and of the next-to-leading order calculation FMNR on the 
right. 
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Figure 9: D*+jet cross sections as function of A<p(D*Jet) compared with the predictions of 
PYTHIA and CASCADE on the left and of the next-to-leading order calculations FMNR and 
ZMVFNS on the right. Due to divergencies in the NLO calculations the last two bins of (a) are 
merged in (b). 
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